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Observed intensities and inferred distances of soft gamma repeaters imply luminosities 
in excess of the nominal (electron-scattering opacity) Eddington limit by four to six orders 
of magnitude. I review the physical basis of this limit. Accretional luminosities may 
exceed it if energy is hydrodynamically coupled from accreting matter to closed field lines 
where it forms a pair gas. This magnetically confined pair gas radiates roughly a black 
body spectrum with ksTe ~ 23 KeV, consistent with observations, at a luminosity up to 
~ 3 X 10'^'^ erg s~^ for a surface field of 10^^ gauss. Magnetic transparency is not required. 
I discuss the minutes-long continuing emission of March 5, 1979, steady counterparts to 
SGR, their spin periods and the recoil problem. 



Subject headings: gamma rays: bursts — stars: neutron — accretion 
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1. Introduction 



The luminosities of the three known soft gamma repeaters (SGR) during outbursts 
are several orders of magnitude greater than the Eddington limit of 2 x 10^* erg of a 
1.4 Mq neutron star with the usual 0.34 cm^ g~^) electron scattering opacity of stellar 
material. For example, the highest luminosity observed from SGR 1806-20 was 1.8 x 10^^ 
erg (Fenimore, Laros and Ulmer 1994), while SGR 0525-66 radiated about 5 x lO^'^ 
erg s~^ (Mazets, et al. 1979) at the peak of its outburst of March 5, 1979. The purpose 
of this paper is to explain how SGR can produce such extraordinarily high luminosities, 
particularly if their energy source is accretion. In §2 I review the derivation and significance 
of the Eddington limit, and discuss its applicability to and implications for three classes 
of SGR models: accretional, thermonuclear and flare. §3 describes how large violations of 
the Eddington limit may be possible in accretional models. In §4 I consider the effects of 
magnetic fields: pair plasma may be magnetically confined by a 10^"^ gauss surface field 
up to luminosities ~ 3 x 10^^ erg s~^, but magnetic transparency is not a satisfactory 
explanation of the super-Eddington luminosities of SGR. Finally, §5 contains a general 
discussion, including an explanation of the extended emission observed on March 5, 1979 
and the implications of the recently discovered steady X-ray sources coincident with two 



Eddington derived his famous limit as a consequence of the equations of hydrostatic 
stellar structure (Chandrasekhar 1939). The equation of hydrostatic equilibrium is 



SGR. 



2. The Eddington Limit 



dp{r) 



GM{r)p{r) 



(1) 



dr 



where M(r) is the mass interior to the radius r, p(r) is the density and p{r) is the pressure. 
An elementary treatment of radiative transfer at large optical depth yields a radiative 
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luminosity 



up dr ' 



where k is the Rosseland mean opacity and Pr = (1 — P)p is the radiation pressure. 
Substituting (1) in (2) gives, for constant (3, the classical Eddington limit Le- 

^ AttcGM ^ AttcGM ^ 

L = (l-/5)< =Le. (3) 

K K 

This derivation begins with the equation of hydrostatic equilibrium; the flow of radi- 
ation is secondary to the contribution of radiation pressure to hydrostatic balance. If k 
were to be increased instantaneously, the consequence would be a reduction in the emer- 
gent luminosity, not the expulsion of matter in a "super-Eddington wind" . Such a wind 
is impossible because there is no regulatory mechanism to hold the radiative luminosity 
constant in an optically thick star. Inefficient convection, an increase of opacity above a 
photosphere, or non-radiative heating of a corona may drive a wind, but these processes 
are outside the scope of Eddington's derivation because the assumptions of hydrostatic 
equilibrium or of radiation diffusion fail; the radiative L does not exceed in optically 
thick regions. 

The Eddington limit is also applicable to accretion flows. Katz (1977, 1987) showed 
that in spherical accretion onto a body of mass M and radius R the rate v^i // of radiation 
diffusion outward through the accreting matter at the surface satisfies 

Vff 6M' 

where the free-fall velocity Vff = {2GM / R)^/"^ and Me = LeR/ {GM) is the accretion rate 
providing the luminosity L^;. A super-Eddington luminosity {L > Le) cannot emerge, even 
if M > Me, because the radiation is trapped by the opaque accreting matter within which 
it is produced. This conclusion was verified in spherically symmetric (Klein, Stockman 
and Chevalier 1980; Burger and Katz 1980, 1983) and disk (Eggum, Coroniti and Katz 
1988) accretion calculations. 
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The Eddington limit is inapplicable to explosions, which are not in hydrostatic equi- 
librium, and which become transparent because their flows diverge spherically. As a result, 
supernovae (L ~ 10'^^ erg s~^) and classical gamma ray bursts (GRB) produced by fire- 
balls at cosmological distances (L ~ 10^^ erg s~-^, the most electromagnetically luminous 
objects in the universe) may exceed their Eddington limits by many orders of magnitude. 
In contrast, the fact that novae have luminosities comparable to Le is evidence that they 
are produced by nearly hydrostatic configurations, rather than by explosions; their mass 
outflow is a consequence of disruption of the erupting star's envelope by its binary com- 
panion, and is incidental rather than essential; if the envelope is not disrupted (because 
the orbit is large) a nova-like variable or symbiotic star is observed. 

No Eddington limit is applicable to energy transported non-radiatively, such as by 
convection in stellar interiors. It is for this reason that the envelopes of cool stars are 
convective; their opacity is very large and the local Le (defined by the local opacity) is 
less than the stellar luminosity. 

An Eddington limit (3), with k > Kes ~ 0.34 cm^ g~^ and Le < 2 x 10^^ erg s~\ 
appears to constrain thermonuclear and accretional models of SGR. The characteristic dy- 
namical time of a neutron star (< 10"'^ s) is much shorter than the observed durations of 
SGR (~ 0.1 s), implying hydrostatic equilibrium (1). Large optical depths (in thermonu- 
clear models because of the required temperature and density, and in accretional models 
because their luminosity requires M ^ Me) imply the diffusion of radiation (2). 

The Eddington limit is not applicable to fiare models, because in them energy trans- 
port is magnetohydrodynamic rather than by diffusion of radiation. The radiating plasma 
must be confined against radiation pressure, but if it is optically thin confinement can be 
provided by a neutron star's magnetic field (Katz 1982, 1993, 1994), without the plasma 
trapping the radiation. Unfortunately, the present understanding of flares is inadequate 
to make specific predictions and fiare models are nearly infinitely fiexible; in addition, 
observed Solar fiares resemble less well the single-peaked profiles and thermal spectra of 
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SGR than the more complex structures and nonthermal spectra of classical GRB. 



In order for an accretion flow to radiate in excess of the accretional luminosity 
must be extracted from the accreted matter itself by a process other than diffusion of 
radiation. I suggest that this may be accomplished magnetohydrodynamically. Suppose a 
discrete body of accreting matter (such as the planetary fragments considered by Colgate 
and Petschek 1981 and Katz, Toole and Unruh 1994) falls onto a magnetic neutron star. 
At impact the accreted matter has speed Vff ^ c/2 and density p ~ 8-100 g cm~^, and 
supplies a dynamic stress pv^^ ~ 10^^-10^^ dyne cm~^. This stress is ~ 3-30 times less 
than the magnetic stress of a 10^^ gauss field, so that if the neutron star has a magnetic 
moment ~ lO^^-lO^-*^ gauss cm^ (as is typical for pulsars) the flow within a few radii of its 
surface is strongly affected by the field. 

The characteristic wave speed in a vacuum magnetosphere, or in one filled with pair 
gas, is ~ c. Because this is only about twice the infall speed, the coupHng efficiency of 
infall kinetic energy to magnetospheric disturbances may be as large as esoR ~ f/c ~ 0.5, 
as would be estimated for an elastic collision between the infalling ffuid and a magneto- 
hydrodynamic wave. Energy coupled to the magnetosphere produces a varying magnetic 
field, and hence an electric field. The amplitude of this electric field may approach that of 
the magnetic field (because the hydrodynamic stress approaches the magnetic stress, and 
the characteristic power density approaches S^c/(87r)), which is large enough to produce a 
vacuum pair breakdown cascade (Epstein and Smith 1993). The detailed electrodynamics 
is complex and beyond the scope of this paper; I simply assume that a fraction esGR of 
the accretional power appears as pair plasma on the closed magnetic field lines, where it 
can be contained by the magnetic stress. 

The Thomson scattering opacity of pair plasma (referred to the pair rest mass only) 



3. Pair Gas and Neutron Star Accretion 



K±= = 

me 




= 726 cm^ g 



-1 



(5) 
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where a^s is the Thomson scattering cross-section, is extremely high, and the corresponding 
Eddington hmit for a 1.4 Mq neutron star would be only 1.0 x 10^^ erg s"-*^. However, the 
particle density of a nondegenerate equilibrium pair plasma with zero net lepton number 
is a sensitive function of temperature: 

n± = p exp(-meC /A;bT). (6) 

The volumetric opacity (beam attenuation per unit length) is 2o"esn±, where the factor of 
two comes from the presence of both e"*" and e~; this is a steeply increasing function of 
temperature for fc^T <^ meC^. 

At high temperatures the particle density of an equilibrium pair gas (6) and its vol- 
umetric opacity are high, and if it is confined the emergent luminosity will be small. At 
lower temperatures n± 0, the plasma is transparent, and its energy density (almost 
entirely radiative) escapes freely even if the charged particles are magnetically confined. 
A photosphere is defined by the usual condition on the optical depth 

2 

T = 2(7esn±X = -, (7) 

where A is a characteristic density scale length. Adopting A = 10^ cm in (7) and using (6) 
yields a characteristic temperature 

Tc « 2.7 X 10^ °K. (8) 

This value of Tc is almost independent of A; dlnTc/dlnX 0.042. 

A nonlinear cooling wave (analogous to an inverse Marshak [1958] heating wave) will 
penetrate a volume of magnetically confined pair plasma, with its visible photosphere at 
an effective temperature Tg ^ Tc- The velocity of this cooling wave is given by energy 
balance arguments as 

VcooliT) = < I, (9) 
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where a is the radiation constant and asB the Stefan-Boltzmann constant. Regions hotter 
than Tc are opaque and have neghgible energy flow; they are hidden from direct observation 
by their own opacity. Cooler regions are transparent; their n± and emissivity are neghgible, 
and they contain chiefly free-streaming radiation. 

The emergent spectrum will be roughly that of a black body at the temperature (8). 
This spectrum is similar to that typically seen in SGR, and explains the observation of 
Fenimore, Laros and Ulmer (1994) that the spectral shape of outbursts of SGR 1806- 
20 is independent of intensity, with only the radiating area varying between outbursts. 
Quantitative predictions of the spectrum require a detailed radiative transfer calculation, 
including the effects of recoil in Compton scattering and cyclotron and cyclotron harmonic 
opacity and polarization, and will depend on the magnitude and direction of the magnetic 
field. The observed peak spectral frequencies of SGR are comparable to the cyclotron 
fundamental or low harmonic frequencies for lO-'^^-lO"'^^ gauss fields, and substantial elliptic 
polarization is likely. 

A characteristic luminosity may be defined: 

Lc = AirR^asBTj ~ 4 x 10^^ ergs"^ > Le- (10) 

SGR need not have luminosities this large; those less luminous may radiate from only 
a fraction of the neutron star's surface, while those with greater luminosity may have an 
effective pair gas photosphere with greater area than the neutron star's surface. The super- 
Eddington luminosity (10) is permitted by the transparency of vacuum; the important 
effect of the magnetic field is to provide a confining stress, preventing rclativistic expansion. 

The harder spectrum observed (Mazets, et al. 1979) during the initial peak of the 
outburst of March 5, 1979 may be explained if during that brief (~ 0.1 s) period the 
pair plasma had too great an energy density to be magnetically confined, and expanded 
relativistically; the transparency wave occurs at the same co-moving temperature T^. (8), 
but the Doppler shift hardens the observed spectrum (Goodman 1986, Paczynski 1986). 
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However, an essential property of the model presented here is that L ^ Le does not 
require an expanding fireball, and only in the first ~ 0.1 s of the burst of March 5, 1979 
was such a fireball likely to have been observed. 



The magnetic field of SGR 0525-66 may be constrained by the observation of a pe- 
riod P = 8 s. If this is interpreted as neutron star rotation, then the radiation rate 
2iJ,'^Q%g/ (3c^) of a rotating dipole fj, (and a similar expression for the wind of the aligned 
component of the dipole moment) lead to the bound on the polar field: 



where t is the neutron star's age, / is its moment of inertia, Ojvs is its spin angular velocity. 



This argument could be applied to the other SGR if spin periods were observed. For 
example, if the suggested (Ulmer, et al. 1993) 2.8 s period of SGR 1806-20 is confirmed 
then (11) would imply Bp < 2.3 x 10^^ (145/^6^4)^/^ gauss. 



Canuto, Lodenquai and Ruderman (1971) calculated the Thomson scattering opacity 
of strongly magnetized plasma with electronic gyrofrequency Og — eB/{meC) in excess 
of the photon frequency uj, and found that for the extraordinary mode the opacity is 
reduced from the unmagnetized value by a factor ~ (a;/f2e)^- A similar reduction factor 
will apply to other electronic contributions to the opacity. The Rosseland mean opacity in 
the presence of a strong magnetic field is therefore ~ 2{u>/Qe)'^ times the zero- field value, 
where the factor of two results from the fact that the magnetic field reduces the opacity 
of only one of the two polarization states. 

Paczyhski (1992) and Duncan and Thompson (1994) pointed out that for sufficiently 
large magnetic fields the observed luminosities of SGR may be sub-Eddington, using the 



4. Magnetic Fields 




(11) 



t4 = t/{10^y), /45 = //(lO^^gcm^), Rq = R/{10^cm), and Pg = ^^/(8s). 



4.1 Magnetic transparency 
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magnetically reduced opacity. The required fields are given by 



LomgC / 2L 



) 



1/2 



B > 



e \Lo 



(12) 



where Lq ~ 2 x 10^^ erg is the zero-field Eddington limit. For photons with hu) = 70 
KeV (near the peak of uFi, for SGR) a luminosity L 1.8 x 10^^ erg s~^ (observed for 
SGR 1806-20; Fenimore, Laros and Ulmer 1994) would require B > 8 x 10^^ gauss, while 
the peak luminosity of 5 x lO"*"* erg s~^ observed from SGR 0525-66 on March 5, 1979 
(Mazets, et al. 1979) would require S > 1.3 x 10-*^^ gauss and its minutes-long continued 
emission with initial L f» 3.6 x 10^^ erg would require B > 1.1 x 10^^ gauss. These 
values are all inconsistent with the bound (11), given the suggested spin periods. 

For very large magnetic suppression of electron scattering opacity it is necessary also 
to consider Thomson scattering by ions. The ionic cross-section 



The corresponding opacity for pure hydrogen khs = 1-16 x 10~^ cm^ g"-*^, while for pure 
iron Kpes — 3.07 x 10~^ cm^ g~^- Assuming planetary cores to be iron, the corresponding 
Eddington limit for a 1.4 Mq neutron star, considering ionic Thomson scattering opacity 
alone (supposing complete magnetic suppression of electron scattering opacity), is 2.3 x 10^^ 
erg s~^, approximately two times too small to permit the observed peak luminosity of 
SGR 0525-66 on March 5, 1979. The ionic opacity may itself be magnetically suppressed 
if fli > a;, but for 70 KeV photons this suppression begins around B 2.4 x lO"*^^ gauss, 
and i? > 4 X 10^^ gauss would be required to explain the observed L ?a 5 x 10'^'^ erg s^-*^. 



The minimum magnetic field required to confine a source of size ~ R radiating a 
luminosity L, continually resupplied, is given by comparison of the radiative and magnetic 
stresses as 




(13) 



4.2 Magnetic confinement 




(14) 
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If the radiation continues over a time r and the radiated energy ~ Lr was not resupphed 
but was all initially contained within a region of volume ~ AnR^/S then the required field 
would be 



For T ^ R/c the condition (14) is less demanding than (15) by a factor which is typically 
~ 100 for SGR. The approximate consistency between the observed duration ~ 0.1 s of 
SGR outbursts and the calculated duration of infall of planetary fragments suggests that 
energy is radiated at the rate at which it is supplied by accretion, so that (14) is applicable. 

Using in (14) black body radiation at an effective temperature equal to the character- 
istic temperature (8) yields a unique value for the required confining field Be- 



a condition easily satisfied by most empirically inferred fields of neutron stars. 

The radiating photosphere may be larger than the neutron star's surface. In this 
case, assuming a dipole for the confining field B, the condition B > Be yields a maximum 
luminosity of confined pair plasma which substantially exceeds that given by (10): 



where B^ is the surface equatorial field. The increase of Lmax over Lc is entirely the 
consequence of the increased radiating area; Tc is unchanged. The peak intensity of SGR 
0525-66 on March 5, 1979 was too large for magnetic confinement unless B^ > 3 x 10^^ 
gauss, as suggested by the interpretation (§3) of its harder spectrum as the result of 
relativistic expansion. All other observed SGR intensities are consistent, by ample margins, 
with magnetic confinement by neutron stars with typical observed pulsar fields. 



The hypothesis of magnetohydrodynamic coupling of accretional energy to a mag- 
netically confined pair plasma can solve the problem of the enormously super-Eddington 




(15) 




(16) 




(17) 



5. Discussion 
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luminosities of SGR. This solution is consistent with their self-absorbed thermal spectra, 
whose temperatures are comparable to that given by (8) and are nearly independent of 
luminosity. Magnetohydrodynamic coupling occurs naturally in accretional models, which 
may be able to explain many of the other properties of SGR. 

At the temperature (8) the thermal velocity of the barycenter of an e+-e~ pair is 
about 0.15c. This suggests the possibility of an observable annihilation line, such as that 
reported by Mazets, et al. (1979). However, they reported it only during the early part of 
the outburst; the relativistic expansion which may have occurred at that time would be 
inconsistent with an observable line. Calculations by Carrigan and Katz (1992) showed only 
a broad bump in the spectrum around 300 KeV, rather than a recognizable annihilation 
line, because in their calculations photons and pairs did not come to thermodynamic 
equilibrium, and the pair temperature was far in excess of that given by (8). They also 
found a nonthermal high energy spectral tail because their assumed source was a power 
law distribution of photons which extended throughout their grid, including transparent 
regions far from the neutron star. 

The 0.1 s intensity peak of March 5, 1979 was followed by a few minutes of continued 
emission with slowly fading luminosity ~ 10^^ erg s~^. This may have been the result of 
continuing accretion, which probably occurred if the distance of closest approach R^. of the 
initial orbit of the accreted fragment to the center of the neutron star lay between R and 
~ 5R. For a one- dimensional distribution of specific angular momentum £ (as expected; 
Katz, Toole and Unruh 1994) the distribution of matter in Rc is dM/ dRc oc Rc dM/ d£ oc 
Rc , while for a two-dimensional distribution dM/dRc is independent of Rc and for a 

1 /2 

three-dimensional distribution dM/dRc oc Rc . In each case the probability of a near-miss 
{R < Rc < 5R) is at least comparable to that of a direct hit {Rc < R). 

The interaction between an accreting disrupted planetary fragment and the neutron 
star's magnetic field is uncertain and complex. In a near miss some matter gives much 
of its angular momentum to the field and is steeply channeled onto the stellar surface, 
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producing the initial spike in intensity. Other matter loses only a little of its kinetic energy 
and angular momentum, and re-emerges from the magnetosphere. Its orbits remain nearly 
parabolic, but their periods are shortened; returning matter arrives in a continuous rain, 
and produces a continuous accretional luminosity. 

If a fraction 6 <^ 1 of the specific angular momentum of an initially parabolic orbit 
is lost by a tangential impulse, the semi-major axis of the new orbit is w i?c/(4(5). The 
new orbital period is 6 x 10~^ {Rc/ {5 R))^/'^ s; a return time of tr corresponds to 5 = 
7 X 10~^ [Re/ R){tr/100s)~'^/^. Because near misses are expected to be about as frequent 
as direct hits, many repetitions of SGR are predicted to be followed by emission of lower 
intensity but longer duration and comparable fluence (Mazets, et al. 1979 estimated that 
in the burst of March 5, 1979 approximately 25% of the total energy was radiated in the 
initial ~ 0.1 s, and 75% in the subsequent ~ 100 s). 

The efl&ciency of conversion of accretional energy to pair plasma is certainly less than 
unity, and very likely less than 0.5. The majority of the accretional energy is trapped in 
optically thick thermal plasma, and will emerge at nearly the Eddington limit. Recent 
observations have found point X-ray source counterparts to SGR 0525-66 (Rothschild, 
Kulkarni and Lingenfelter 1994) and SGR 1806-20 (Murakami, et al. 1994; Sonobe, et 
al. 1994) with steady luminosities 10^^ erg s~^ and f=i 3 x 10^^ erg s~^ respectively. 
These steady luminosities exceed the expected luminosities of cooling neutron stars of the 
estimated ages, and it is possible that they are associated with the accretional energy 
released in the outbursts. However, they are well below the nominal (electron scattering) 
Eddington limit, and thus apparently inconsistent with their interpretation as the gradual 
release of trapped accretional energy. 

The steady luminosity of SGR 1806-20 exceeds its time-averaged burst luminosity ~ 
100-fold. The steady luminosity of SGR 0525-66 exceeds its time-averaged burst luminosity 
(excluding March 5, 1979) by an even larger factor but is comparable to the energy released 
on March 5, 1979 divided by the time since that date. These facts also suggest that the 
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trapped accretional energy may be insufficient to power the steady counterparts. 

The classical power source of neutron stars is their rotational energy. If the apparent 
8 second period of SGR 0525-66 is its spin period, then its spin-down power IVl\ig/{2t) ~ 
10^^ /t4 erg s~-^, inadequate to explain its observed ~ 10^^ erg s~-^ steady counterpart. A 

— 1/2 

spin-down power supply would be adequate if P < 0.3^4 s, or if the SNR contains a 
second neutron star with the required smaller P. The steady X-ray counterpart to SGR 
1806-20 of ~ 3 X 10^^ erg s~^ similarly would require P < 0.5t^^^^ s, consistent with all 
known data (except the very uncertain suggestion of a 2.8 s period by Ulmer, et al. 1993). 

The requirement that the present spin-down power be supplied by a neutron star 
of age t may be used to constrain its magnetic field. The neutron star's spin rate is 
eliminated from (11) by requiring that the observed X-ray luminosity Lx = lO^^Lae erg 
s~^ be powered by spin-down with efficiency ex = O.Oleo.oi- This yields a bound on the 
magnetic dipole moment 



The presence of neutron stars with comparatively fast spin {P < 0.3 s) would also 
explain the observation that the SNR are "plerionic" (center-filled) and apparently powered 
by a neutron star rather than by the collision of stellar debris with interstellar gas. The 
methods of estimation of age applicable to collision-powered SNR (Kulkarni, et al. 1994; 
Rothschild, Kulkarni and Lingenfelter 1994) would then be inapplicable, and their actual 
ages could be significantly longer. If t > 10^ y then the inferred neutron star recoil 
velocities would be consistent with the retention of planets, as required in models of SGR 
which assume accretion of planetary fragments, even if the SGR are dynamically off-center 
in their SNR. 

I thank NASA NAGW-2918 for support and N. G. H. Sh. for illuminating these 
apparently miraculous events. 
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